Problems With Extracting m s from Flavor Breaking in Hadronic r Decays 
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A numerical error is pointed out in the existing expression for the 0(o? s ) longitudinal component 
of the squared-mass (D — 2) contribution to the hadronic r decay rate. The corrected version is 
found to be such that, to O(af), each term in the resulting series is larger than the previous one, 
hence ruling out the direct use of flavor breaking in hadronic r decays as a means of extracting 
the strange quark mass. An alternate approach, in which one uses the model spectral functions 
previously developed for the strangeness-changing scalar channel (and employed in alternate sum 
rule analyses of m s ) as input to the r decay analysis, is shown to provide mutual consistency checks 
for the two different methods. The results for m s implied by r decay data, and the existing model 
spectral function are also presented. 
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As is well-known, the inclusive hadronic r decay ratio 

r[r~ — > v T hadrons (7)] 



R T = 



T[t -> v T e v e {l)} 



(1) 



(where (7) indicates additional photons or lepton pairs) may be rather reliably computed using techniques based on 
the OPE and perturbative QCD (pQCD) [1-10]. (A very clear recent review of the theoretical situation is given in 
Ref. [11].) One begins with the representation of this ratio in terms of hadronic spectral functions, 



R T = \1itSew 



ds 



1 + 2-^- J ImnW(s)+Imn(°)(s) 



(2) 



where Sew — 1-0194 represents the leading electroweak corrections [12], J = 0, 1 labels the hadronic rest-frame 
angular momentum, and the W- J '(c[ 2 ) are given in terms of the vector and axial vector current spectral functions 



n^v^O? 2 ) ( with y = ud,us) by 



n«( s ) = \v ud \ 2 K J )( S ) + n« ( S )l + \v us \ 



n 



us,V 



(3) 



where 



i J ^ e i ^{0|T(J4 V)j4 (a ; )J-,^(0)t)|0) =(-^ g 2 + ^)ng^ A (g 2 ) + QV< ; V,A(9 2 ) 



(4) 



with J^j.y A the standard vector and axial currents involving flavors ij, defines H^'p v A (q 2 ). Using the analyticity 
properties of the correlators, the hadronic representation, Eq. (2), can then be converted into a contour integral 
representation 



R T = 6ttSew i <t> —5 ( 1 ~o 



,2 mt 



i + 2^)n( 0+1 )( s )-2^n(°)( s ) 



m 2 . J 



(5) 
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Since m 2 ~ 3 GeV 2 , the latter expression is amenable to evaluation using the OPE. Keeping operators up to dimension 
D = 6 in the OPE, and evaluating the mass-independent (D = 0) and mass-dependent (D = 2) pcrturbative terms 
to 0(o? s ) and 0(a 2 s ), respectively, it is known that an excellent fit to hadronic r decay data is obtained using a value 
of a s (m T ) compatible (after running) with that obtained directly at a scale Mz (see Ref. [11], and earlier references 
cited therein). Of relevance to the present paper is the observation, made by numerous earlier authors, that, although 
the D = 2 ud contributions to R T are tiny, the D = 2 us terms can alter the prediction for the us contribution to 
R T by ~ 20%, for conventional values of m s . (See, for example, Ref. [4,6,13,14].) These contributions, of course, play 
little role in determining R T , owing to the suppression by the factor of \V US /V u d\ 2 relative to the non-strange current 
contributions. As noted by Davier [13], however, sufficiently precise r decay data allow one to look explicitly for this 
effect, and hence, in principle, determine m s . 

The most convenient way of implementing this extraction of m s is to consider the contributions to inclusive hadronic 
r decay mediated separately by the strangeness-changing (us), and strangeness-non-changing (ud) weak currents, and 
take the difference of these contributions after rescaling each by the inverse of the square of the corresponding CKM 
matrix element. The D = contributions then cancel in the difference, leaving a leading D = 2 contribution 
proportional (to leading order in m u /m s ) to m 2 . Numerical estimates show that this term is not only formally 
leading, but also numerically leading at the scale m T fixed by the radius of the circular contour in Eq. (5) [4,6]. 

The extraction of m s , as performed by the ALEPH Collaboration [13,14], using flavor breaking in hadronic r decays, 
relies on the expressions for the D = 2 mass-dependent terms worked out by Chetyrkin and Kwiatkowski [6] . Having 
determined the D = 2 contributions to the correlators n^^^ ^(g 2 ) and performed the contour integration indicated 
in Eq. (5), they obtain the result, to leading order in m u /m si quoted in their Eq. (26). This result implies, summing 
"longitudinal" (J — 0) and "transverse" ( J = + J = 1) contributions, a us current contribution to R T given by 



i£>=2 

t\ U s-,v,a 



T2 \V U , 



~EW ' 



m 2 (m T ) 



mt 



16 

1 + — a(m T ) + 11.03 a(m T ) 2 



(6) 



In Eq. (6), m s (m T ) is the running strange quark mass at a scale m T , and a(ji) = a s (fi)/n, with a s (n) the standard 
MS running coupling at scale /i. The polynomial P(a) = 1 + 16a/3 + 11.03a 2 appearing in Eq. (5) is evidently rather 
well-converged to 0(a 2 ), allowing, once experimental data is sufficiently precise, a reliable determination of m s (m T ). 
However, as we will see below, the coefficient 11.03 multiplying a 2 in P(a) is incorrect. Not only is the corrected 
value (46.002) such that the convergence of P(a) is not very compelling (P(a) = 1 + .602 + .587+ ... at the scale m T ) 
but, worse, this large value is generated mostly by an enormous 0(a 2 ) coefficient in the longitudinal contribution. 
The series for the longitudinal contribution then turns out to be such that each subsequent term is larger than the 
previous, making it impossible to rely on the OPE representation of this contribution, and hence impossible to use 
the OPE representation of the D=2 terms to extract m s . 

We will, however, see that it is still possible to make some use of the data on flavor-breaking in r decay, since 
the spectral function of the badly behaved longitudinal contribution is related to that occuring in existing sum 
rule analyses of m s based on the strangeness-changing scalar channel. Although the spectral function has not been 
experimentally determined in that channel, model versions exist [15-17], and information on the low-energy portion of 
it is available using the value of the scalar Kn form factor measured in together with the Omnes representation 
of the form factor, which takes experimentally measured Kir phase shifts as input [17]. 

The rest of this note is organized as follows. In Section II we present the corrected results for the longitudinal and 
transverse D = 2 us contributions to R T and demonstrate the non-convergence of the longitudinal series in a(m T ). In 
Section III we will see how to reformulate the expression for the us contribution to R T in such a way that it involves 
only (1) a longitudinal contribution which can be calculated directly from any given (model-dependent or otherwise) 
version of the spectral function for the strangeness-changing scalar channel, (2) reasonably well-converged D = 2 
transverse contributions and (3) higher dimension longitudinal and transverse contributions which are either quite 
accurately known or small. In Section IV we will then use existing models for the scalar spectral function, together 
with recent r decay data, to perform a consistency check on the values of m s extracted using the two different analyses. 



II. THE PERTURB ATIVE SERIES FOR THE D=2 CONTRIBUTION 



The D = 2 us contributions to R T follow immediately from Eq. (5) and the expressions for the longitudinal and 
transverse pieces of the relevant vector and axial vector current correlators to next-to-next-to-leading order, as given 
by Chetyrkin and Kwiatkowski [6]. In what follows we will drop terms of order m u /m s since, as we will see below, it is 
possible to employ r decay data (once we are aware of the poor convergence of the longitudinal contributions) only if 
we neglect such contributions. Note also that the same extremely poor convergence of the longitudinal contributions 
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also afflicts the D = 2 ud contributions; however, in this case, the squared-mass factors are so tiny that this is of no 
numerical significance in existing analyses of the isovector hadronic weak current mediated decays. 
The contributions to the correlators, in this approximation, are then given by [6] 



n (0)D=2 
ll us;V,A 



3 m s (Q 2 
2tt 2 Q 2 



+ 



(7) 



and 



n (0+l)D=2 



3 m s (Q 2 ) 
4tt Q 2 



1 + -a(Q 2 ) + 19.5831a(Q 2 ) 2 



(8) 



where a(Q 2 ) = a s (Q 2 )/ir and m s (Q 2 ) are the running coupling and running strange quark mass, both at scale 
/i 2 = Q 2 = — s, in the MS scheme. In Eq. (7), + • • • stands for the second set of terms of Eq. (16) of Ref. [6], which 
terms do not survive the contour integration involved in obtaining R T . Note that, in Eq. (8), the coefficient 19.5831 
includes the contribution of the three-loop graph containing two quark loops, which contribution is also present in 
the corresponding ud contribution to R T and hence cancels in forming the flavor-breaking difference discussed below. 
(The term which survives corresponds to the replacement 19.5831 — > 19.9332.) Note that there is a typographical 
error in Eq. (25) of Ref. [6], where the correct result 19.5831 ~ 19.6 has inadvertantly been recorded as 17.6. 

It is now a straightforward matter to evaluate the D = 2 us contributions to R T . One first expands a(Q 2 ) in terms 
of a(m 2 ) = a and m s (Q 2 ) in terms of a and m s (m 2 ) = m s , the coefficients in these expansions being polynomials in 
log (Q 2 /to 2 ) [18]. The remaining integrals are then elementary ones, involving products of powers of Q 2 and powers 
of log (Q 2 /m 2 ). Performing these integrals, one finds, for the contributions to R T corresponding to the terms written 
explicitly in Eqs. (7) and (8) above, 



4 0) 



D=2 



— —3\V US \ Sew 



m 2 



us;V,A 

D=2 _ 



1 + —a+ 109.9889a 2 



us;V,A 

whose sum yields, as claimed above, 



= -9|K a | 2 S w -f [1 + 4a + 24.6707a 2 ] , 



rn 



i R r]us:V,A - ~^\V us \ 2 S E w 



^2 r 



16 , 
1 + — a + 46.002a' 



(9) 
(10) 

(11) 



(In the transverse contribution, 24.6707 — > 25.0207 when one takes into account the cancellation of the terms propor- 
tional to a 2 m 2 common to the ud and us correlators, and associated with the graphs involving a quark bubble on the 
internal gluon line.) The revised numerical values quoted above have been confirmed by Chetyrkin [19]. 

If one uses the value of a s (m T ) determined from the ALEPH analysis of non-strange hadronic r decays [20] (this 
choice being the most sensible since, in contrast to the global analysis, the results are not affected numerically by the 
use of the erroneous form of the D = 2 contributions), one finds, for the polynomial appearing in the Eq. (9), 



1 + -8a + 109.9889a 2 



1 + 1.052+ 1.397 . 



(12) 



It is quite clear that the series is simply too poorly converged to allow one to consider the coefficient of m 2 appearing 
in the expression for R T reliably determined by the OPE analysis. As such, the conventional approach to obtaining 
m s from flavor breaking in hadronic r decay is ruled out. 



III. A MODIFIED ANALYSIS FOR m s 

Two observations make it possible to utilize flavor breaking in hadronic r decay to constrain m s , in spite of the non- 
convergence of the series for the longitudinal contribution noted above. The first is that, to the extent that one ignores 
corrections of order m u /m s , the vector and axial vector D = 2 us contributions are identical. (Beyond leading order 
in m s , this statement is, however, no longer valid.) One may then, as far as the D = 2 contributions are concerned, 
concentrate on, say, the vector contributions. This is of relevance because of the second observation, namely that the 
vector current spectral function is related to the spectral function of the strangeness-changing scalar channel, about 
which some (though not complete) experimental information is available. Combining these two observations, one sees 
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that the problematic D — 2 contributions can be handled by replacing the full (D = 2,4,6, ... and vector plus axial 
vector) longitudinal contributions with twice the full vector contribution, providing one makes corrections reflecting 
the difference of the vector and axial vector contributions. These corrections then (to leading order in m s ) begin at 
D = 4. The full vector contribution is obtained, as described below, from the model spectral functions employed in 
the strangeness-changing scalar channel. 

Let us define R\^ v °a to be the vector (V) or axial vector (A) flavor ij = ud or us contributions to R T associated 
with those operators of dimension D in the OPE of the longitudinal (J = 0) or transverse (J = + 1) terms in 
Eq. (5) after scaling out the factors |Vy | 2 Sew, an d Rij to be the corresponding sum, for given ij, over D, V + A 
and J 



0,0 + 1. As noted above, to be able to effectively employ the model spectral functions from the scalar 
channel as input it is necessary to drop corrections of 0{m u /m s ) to R^I-va 2 - F° r consistency one must then also 
drop other terms proportional to the light quark u, d masses. For D = 2, this includes all of R^lyA 1 apart from 
the "quark bubble" transverse contribution proportional to a 2 m 2 mentioned above, which simply cancels against the 
corresponding term in R^'y^ 2 - F° r D = A, since the leading contributions are proportional to < m s ll >, I = u,d, 
and m,g, this means dropping also terms proportional to < mfll >, and terms fourth order in the quark masses, 

other than those proportional to mf. This removes all contributions to R^d-VA^ a P ar t from those terms in the 
transverse contribution proportional to the gluon condensate, which terms cancel in forming the difference of ud and 
us contributions. Note also that, to the order so far computed, the D = 6 contributions are pure transverse. With 
the set of approximations just discussed, one thus obtains, 
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Ry 
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(0);D=4 



R 



(0);D=4 



R 



(0+l);D=2 

<;V 



+ R 



(0+l);£>=2 
us; A 



x \ " r R (°+ 1 );- D , p (0+i);O _ p(o+i);£> _ p(o+i);D 

[-"-ud-.V ' n ud;A n us;V n us;A 



D=4,6,- 



(13) 



where terms with D > 6 have been dropped. The first term on the RHS of Eq. (13) is to be obtained using the 
input model spectral function obtained from the scalar channel. The remaining terms will be obtained from the OPE 
representation. The prime on the third term in Eq. (13) reminds us that the "strange quark bubble" contribution 
cancelled by the corresponding contribution to the ud correlator has to be removed in evaluating this term. The ud 
contributions in the last term serve only to remove all D = 4 contributions proportional to the gluon condensate. As 
we will see below, all of the remaining terms are either very well, or relatively well, known. Eq. (13) thus isolates 
the problems of the D = 2 part of the longitudinal term in the OPE representation in such a way that, to the extent 
one can obtain a good phenomenological representation of the scalar spectral function, one can avoid the poorly 
converging OPE representation and instead evaluate this contribution phenomcnologically. While, at present, it is 
not possible to do this in a completely satisfactory manner (see the discussion on the existing constraints on the scalar 
spectral function below), the situation is amenable to future improvement. 

Let us turn then to the evaluation of the remaining terms in Eq. (13). We stress again, that in writing down all the 
expressions below, corrections proportional to the light quark masses have, for consistency's sake, been systematically 
dropped, with the exception of one term for which there is a numerical enhancement, as discussed below. Since the 
expressions for the relevant contributions to the correlators are well-known in the literature [4,6] , we present only the 
OPE representations after the relevant contour integrations have been performed. We then find, 



R (0);D=4 _ „(0);D=4 
n us;A n us;V 



24tt 2 



[-2m s <£l>]+ 12 



12 5 
Ta + 7 



rris 



mi 
m s 



R (0+1);D=2 p (0+l);D=2 



' = -18 









[1 + 4a + 25.0207a 2 + / 3 a 3 ] 



R 



(0+l);D=4 + ^(0+l);D=4 _ p (0+l);D=4 _ D (0+1);D=4 



ud-V 



ud;A 



-R 



us\V 



-R 



-54 



m-r 



us\A 
4 



R 



(0+l);U=6 + R (0+1);D=6 _ R (0+1);D=6 _ R (0+1);D=6 



ud-y 



ud:A 



R 



'us;V 



R 



us:A 



81^ 2 a 2 . n 

2 — -2m s < ss > 



. 



(14) 



A few comments are in order concerning Eqs. (14). First, the leading term fourth order in the quark masses in the 
first expression, though formally an 0{mi/m s ) correction to the contribution appearing in the third expression, 
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has a numerical enhancement due to the presence of a factor 1/a in the coefficient multiplying it, which enhancement 
largely undoes the mi/m s suppression. This term has, therefore, been kept, though it turns out that, in fact, 
both terms are numerically tiny, and play a negligible role in the final analysis. The presence of a 1/a factor in 
the coefficient multiplying m 3 mi is due to the well-known fact that the quark condensate terms obtained from the 
standard Wick's theorem reduction must be modified to absorb additional long-distance mass logarithm terms. The 
modified condensates are then no longer scale invariant. Defining new scale invariant versions (see Rcf. [4] for details) 
one finds that, rewriting the D = 4 contributions in terms of these condensates, such inverse power dependence on 
a shows up in the modified fourth order mass terms. The quark condensates appearing in Eq. (14) are thus those 
defined in the Appendix of Refs. [22,4]. Second, in the D = 2 transverse terms of the second line, the coefficient 
f 3 of a 3 has not, as yet, been computed. Because of the slow convergence of the 0(a,a 2 ) terms in the series, we 
have, however, made an estimate of / 3 . This is done by first estimating the corresponding coefficient in the scries 
for the D = 2 contribution to the correlator using the method of Ref. [23], and then integrating the result around 
the circular contour. The proceedure of Ref. [23] is known to produce values for the 0(a 3 ) coefficient accurate to 
±25 in the case of those mass-dependent observables for which the 0(a 3 ) terms have been computed [23]. We find 
fs = 217.8 from the "fastest apparent convergence" version of the method, and = 219.6 for the "principle of minimal 
sensitivity version". As a conservative estimate of the effect of higher order terms we thus take fz — 220 ± 220. At 
present the uncertainties in fh s that result from that in fa are much smaller than those associated with the errors 
in the experimental input. Finally, the reason for setting the D = 6 (purely transverse) contributions to zero needs 
explanation. It is straightforward to read off the explicit forms of these contributions from Eq. (3.12) of Ref. [4]. As 
discussed in that reference, however, the full set of D = 6 condensates is not known empirically, and it is customary 
to make an estimate based on a rescaled version of the vacuum saturation hypothesis (Eq. (3.13) of Ref. [4]). This 
estimate has relatively large errors, but can be checked by employing the method of spectral moments [5,11], which 
allows one to make an empirical extraction of (at least the ud) D = 6 contributions using the measured spectral data. 
This extraction has been performed by the ALEPH collaboration [20,21], who obtain, for the sum of vector and axial 
vector ud D = 6 contributions, a value seven times smaller than the central value obtained in Ref. [4], but compatible 
with it, within the estimated theoretical errors. With this empirical information as input, and the realization that 
the difference of the ud and us D — 6 contributions should be further suppressed by the approximate SU(3)f flavor 
symmetry, the numerical value of the combined D = 6 contribution, though not well-determined, becomes completely 
negligible numerically, and hence is neglected in the analysis which follows. 



IV. A MODIFIED EXTRACTION OF m. 



In what follows we detail the input required to numerically evaluate the contributions listed in the previous section. 
We employ, from the various ALEPH fits for a = a s (m T )/TT, that obtained from the analysis of the non-strange 
modes alone, since the slow convergence of the longitudinal contributions in the strange case make that extraction 
unreliable. In quoting values for m s , we will also run the value extracted at scale m T down to a scale 1 GeV 2 , using 
the recently-determined four-loop (3 [24] and 7 [25] functions, in order to provide comparisons to other results for 
m s [15-17,26,27] all of which are quoted at that lower scale. Other necessary inputs, apart from the longitudinal 
vector contribution, to be discussed in detail below, are the light quark mass ratio m s /m = 24.4 ±1.5 [28] (where 
to = (m u + TO,j)/2) and 

- 2m s < ss >= (^) ( <J*> ) (-2m < uu » = tt) ( <Ji> ) /> 2 , (15) 
\m J \<uu>J \m ) \<uu>) n n 

where, following Refs. [15,16], we take the ratio of strange to light quark condensates to lie between 0.7 and 1. 

The crucial piece of input, however, which can, in principle at least, allow us to evade the problems with the longitu- 
dinal D = 2 terms, is the full longitudinal vector contribution, ^2 D R^s-'y • This can be obtained straightforwardly 

from any model (or empirical determination) of the spectral function of the strangeness-changing scalar channel since 
the longitudinal us vector current spectral function is simply 1/s 2 times that of the scalar channel. The latter spectral 
function begins at Kir threshold, the next open channel being Kiririr (experimentally, however, the Kir channel is 
purely elastic below (1.3 GeV) 2 [29]). Given that the K*(U30) has a branching ratio of 93 ± 4 ± 9% to Kir [30], it 
is safe to assume that the Kir mode dominates the scalar spectral function out to s ~ 2 GeV 2 . This observation is 
particularly relevant because the Kir contribution to the spectral function is determined by the timelike Kir scalar 
form factor, which satisfies an Omnes representation [15]. Using as input the experimental value of the form factor 
determined in Kg 3 , together with the measured Kir phase shifts, it is thus, in principle, possible to determine the 
Kir contribution to the physical scalar spectral function [15,17]. In practice, since the Kir phases have only been 
measured out to s = (1.7 GeV) 2 , certain assumptions are required about the behavior of the phase beyond this point. 
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In Ref. [17] it has been assumed that the asymptotic value, ir, of the phase (associated with the known asymptotic 
1/s-dependence of the scalar form factor) is achieved for all s > (1.7 GeV) 2 . At present, neither theoretical nor 
experimental checks of this assumption exist, though the experimental phase does approach the asymptotic value as 
s — ► (1.7 GeV) 2 . One may also check that, if for example the phase were to make an excursion into the third quadrant 
as one passed through the Xq(1950) region, and then return to ir, the longitudinal spectral integral relevant to r 
decay would be increased by less than <~ 10 — 15%, producing a decrease of m s of less than ~ 2%. Were the onset of 
the asymptotic behavior to occur significantly beyond this point, however, the scalar form factor obtained in Ref. [17] 
could be significantly altered, even at low s. This difficulty is, of course, potentially surmountable in future. Lacking 
any additional information at present, and in light of the discussion above, we will tentatively accept the assumption 
of Ref. [17], and hence consider the Kir portion of the spectral function, and hence the full scalar spectral function 
below ~ 2 GeV 2 , to be as determined there. 

Even if this assumption is correct, however, this leaves the problem that one expects additional channels to become 
increasingly important as s is increased beyond the location of the Kq(1A30). Since the ^(1950) Kir branching 
fraction is only 52 ± 14% [30], it is likely that employing only the Kir contribution to the full spectral function 
will become a questionable approximation significantly before one has reached the end of the range of the spectral 
integral relevant to r decay. This likelihood is reinforced by the observation that, attempting to perform a standard 
Borcl transformed QCD sum rule analysis of the scalar channel using only the Kir portion of the spectral function 
determined in Ref. [17], one finds no stability plateau for the "extracted" strange quark mass in the range s ~ 3 GeV 2 
[31]. Since 




where Pu S ( s ) ls t ne spectral function of the scalar strange channel, such an underestimate of the spectral function 
corresponds to an underestimate of the full longitudinal contribution to r decay, and hence (because of the overall sign 
(+) with which the contribution in Eq. (16) enters R uc i — R U s), to an overestimate of m s . While the weight function 
entering the spectral integral in Eq. (16) has a double zero at s = m 2 , this does not help the situation as much as 
one might hope: although the maximum of the weight function occurs for s/m 2 = 1/3, the weight has decreased by 
only a factor of 2 from its maximum value for s/m 2 = 2/3. At present there is, thus, no satisfactory extension of 
the full spectral function beyond s ~ 2 GeV 2 , and this limits our ability to extract m s from the hadronic r decay 
data. Note, however, that, while neglect of the non-i^Tr portion of the spectral function leads to an overestimate of 
m s when one employs the r decay analysis, it leads to an underestimate of m s when one employs the direct QCD sum 
rule analysis of the scalar channel. The extent to which the two different extractions of m s (both of which depend on 
the longitudinal spectral function) agree thus provides some post facto information of the degree of reliability of the 
model spectral function. 

Using the Kir portion of the spectral function as determined in Ref. [17], we then find 

= .055 . (17) 

D 

One can only guess at the error on this number at present. Varying the parameters of the fit to the experimental 
phases as described in Ref. [17] one finds a variation of ±.0012, while, as explained above, the uncertainty associated 
with the unknown behavior of the Kir phase above (1.7 GeV) 2 could easily produce uncertainties of order 10 — 15% 
or more. By way of contrast, the model spectral function employed in Refs. [15,16] produces a value .090. Note, 
however, that, as pointed out in Ref. [17], the behavior of this spectral function below s ~ 2 GeV 2 is incompatible 
with the Omnes representation and known Kir phases, and hence is unphysical. 

For the experimental input to the analysis we employ the following. First, for the strange and non-strange current 
contributions to R T in Eq. (1), we employ the latest ALEPH results [13,14,21,34] 

R T . ud = 3.493 ± .026 

R T . US = 0.155 ± .008 . (18) 

These are compatible with those of CLEO [32,33]; the latter collaboration, however, has not yet released an official 
number for the strange branching fraction [33]. For the CKM matrix elements, we take the value of V us from K e3 
data, where the extraction is under best theoretical control, and V u d from three-family unitarity (there exist, for 
example, effective isospin-breaking electromagnetic contact interactions, which involve no explicit photons in the 
effective hadronic theory, and which are, therefore, not taken into account in treatments of super-allowed (3 decays; 
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it is not clear how to determine a probable error associated with this neglect). With this input, and the electroweak 
correction from Ref. [12], we obtain 

[Rud - Rus] expt = 0.446 ± 0.169 . (19) 

The cancellation present in Eq. (19), of course, magnifies the errors. For example, if we employed the older ALEPH 
value R T -us = 0.156 the central value would be shifted from 0.446 to 0.426. The overall error quoted is dominated by 
that on R us . 

From Eqs. (13), (14), (15), (17) and (19), one obtains the following upper bound for m s 

m s (l GeV 2 ) < 220t!$ ± 33 MeV » ( 20 ) 

where the first set of (asymmetric) errors is associated with the experimental uncertainties in Eq. (19), and the second 
with the estimate above of the uncertainties due to higher order contributions in the perturbative series for the D = 2 
transverse terms. 

In quoting errors in Eq.(20) we have refrained from including those due to possible shortcomings in the model 
longitudinal spectral function (associated with (1) the unknown high energy behavior of the Kir phase and (2) 
neglected Kmrir, ■ ■ ■ contributions to the spectral function). We remind the reader that these corrections, if made, 
would produce a value lower than the bound shown in Eq. (20), though it is not possible to make a sensible estimate 
of the size of this correction at present. In contrast, the same set of corrections would serve to raise the value 

125 MeV < m s (l GeV 2 ) < 160 MeV, (21) 

obtained from the analysis of the scalar sum rule in Ref. [17]. The true result for m s will, in general, lie between the 
values obtained from the r decay analysis and the scalar sum rule, so long as one uses only the Kir portion of the 
longitudinal spectral function in both analyses. At present, the two results, owing to the large errors, especially in 
Eq. (20), are compatible, but one would like to see the errors in both reduced. 

To conclude, we stress that flavor breaking in hadronic r decays can still be employed as a means of extracting 
m s , despite the poor convergence of the perturbative series for the D = 2 longitudinal terms, provided empirical 
input can be obtained for the longitudinal spectral function. As pointed out in Rcfs. [35,36], it is possible to separate 
the longitudinal and transverse components experimentally by analyzing the dependence of the cross-section on the 
angle, /3, between the direction of the K and the laboratory, as seen from the hadronic rest frame. One can thus, 
using t decay data alone, obtain an estimate for m s , in the approximation that one neglects the higher multiplicity 
(Kmnr, ■ ■ •) contributions to the longitudinal spectral function. Since measuring the Kir component of the longitudinal 
spectral function would simultaneously determine the Kir component of the spectral function for the scalar channel, 
a comparison of the m s values extracted using r decay and the scalar channel sum rule would set upper and lower 
bounds for m s , and provide an explicit check of the reliability of the the approximation of neglecting higher multiplicity 
intermediate states. The experimental determination of the Kir component of the longitudinal spectral function is, 
of course, also of interest in that it would provide implicit information on the high energy behavior of the Ktt phase. 
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